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摘要 :目的 提出 一 种 结合 低 秩 与 稀疏 惩罚 的 PET 动 态 图 像 有 


方法 (L&S)。 方 法 建立 L&S 重 建 模 型 ,利用 split Bregman 法 


来 最 优化 求解 代价 函数 。 采 用 单 房 室 模 型 仿真 一 套 PET 心 肌 82Rb 港 注 图 像 , 将 L&S 重 建 方法 与 最 大 似 然 期 望 值 法 (MLEM)、 
低 秩 惩罚 和 稀 焉 惩罚 重建 方法 比较 。 结 果 L&S 方 法 重建 的 图 像 的 均 方 误差 (MSE) 最 小 ,并 且 保 留 了 更 多 图 像 特 征 。 男 外 


L&S 重 建 得 到 的 靶 心 图 和 参考 组 的 靶 心 图 最 相近 。 结 论 LSH 


J 
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让 建 方法 无 论 是 在 直观 视觉 上 ,还 是 定量 分 析 上 都 优 于 男 外 3 种 
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Abstract: Objective To propose a new method for dynamic positron emission tomographic (PET) image reconstruction using 


low rank and sparse penalty (L&S). Methods The L&S reconstruction model was established and the split Bregman method 
was used to solve the optimal cost function. The one-tissue compartment model was used to simulate a set of PET 82Rb 
myocardial perfusion image. The L&S reconstruction method was compared with maximum likelihood expectation 


maximization (MLEM) method, low-rank penalty method and sparse penalty method. Results The L&S reconstruction 
method had the smallest MSE and well maintained the feature information. The polar map created by L&S method was the 


most similar with the reference actual polar map. Conclusion L&S reconstruction method is better than the other three 
methods in both visual and quantitative analysis of the PET images. 
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与 传统 的 SPECT 心 肌 成 像 相 比 ,PET 心 肌 灌注 成 
像 可 以 提高 诊断 准确 性 ,更 好 地 确认 病灶 位 置 ,还 能 定 
量 分 析 心 脏 血 流 沁 。 尤 其 是 PET 动 态 心肌 灌注 成 像 还 
能 测量 心肌 区 域 放射 性 药物 的 生物 分 布 呈 ,这 就 使 得 估 
计 动 力学 参数 [例如 示 踪 剂 输送 速率 参数 K 和 心肌 血 
流量 (MBF) ] 变 得 非常 有 意义 。 然 而 ,PET 动 态 心肌 灌 
注 成 像 主要 还 处 于 研究 阶段 ,临床 上 还 缺乏 应 用 。 主 要 
的 问题 是 将 数据 细 分 到 短 帧 时 ,会 华人 噪声 ,从 而 影响 
心肌 绝对 血 流 的 定量 。 传 统 上 ,标准 的 动态 PET 成 像 由 
单独 的 帧 独立 重建 得 到 ,然后 在 体 素 或 ROI 水 平 上 运用 
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动力 学 模型 得 到 时 间 活 度 曲 线 (TAC)。 独 立 图 像 重 建 
主要 靠 统 计 图 像 重 建 方法 实现 ,例如 最 大 似 然 期 望 值 法 
(MLEM) 中 。 然 而 , 当 在 低 计数 情况 下 ,直接 采用 
MLEM 算 法 方差 很 大 "。 而 且 方 差 还 会 随 着 采样 间隔 
的 增 大 而 增 大 。 有 人 提出 利用 PET PERS i PE HY 
重建 算法 ”。 例 如 在 PET 动 态 图 像 重建 中 加 入 全 变 分 
正则 约束 。 这 种 方法 不 论 是 在 图 像 空间 哺 , 还 是 在 测度 
空间 所 上 都 得 到 了 应 用 。 另 外 多 维 小 波 去 噪 可 以 恢复 
隐藏 在 动态 PET 图 像 中 的 生物 信号 保 真 度 ,从 而 准确 的 
定量 心肌 灌注 。 学 者 Su 和 Shoghi 在 2008 提 出 小 波 降 
噪 技术 ", 它 对 噪声 敏感 度 低 ,在 高 噪声 水 平 下 可 以 提 
供 更 加 准确 的 参数 估计 。 

最 近 , 有 研究 者 提出 了 利用 矩阵 的 低 秩 性 质 恢复 图 
像 。 该 研究 首先 在 动态 核磁 共振 成 像 汪 ”上 受到 关注 ， 
学 者 Lingala 等 “将 其 与 稀 玻 先 验 结合 分 别 运用 到 磁 
共振 图 像 的 恢复 和 重建 上 。 在 动态 PET 心 肌 灌注 成 像 
上 ,连续 的 图 像 之 间 有 高 度 相 关 的 元 余 信 息 , 这 使 得 网 
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像 矩 阵 有 很 好 的 低 秩 性 。 另 外 在 图 像 某 些 部 分 ,含有 血 
流 灌注 信息 ,这 部 分 是 稀 踊 的 。 因 此 ,本 文 在 动态 PET 
心肌 灌注 成 像 上 提出 一 种 结合 低 秩 与 稀 玖 惩罚 的 重建 
方法 (L&S)。 我 们 利用 split Bregman 算 法 求解 凸 最 优 
化 问题 ,并 利用 真实 的 PET 心 肌 容 b 灌 注 成 像 来 最 优化 
正则 化 参数 。 我 们 在 成 像 视觉 上 和 信 噪 比 上 对 L&S 方 
法 做 了 评估 。 同 时 我 们 将 该 算法 与 传统 的 MLEM 算 法 
以 及 单独 利用 稀 玻 或 低 秩 惩 罚 的 算法 进行 了 比较 。 


1 材料 和 方法 
1.1 动态 2RDPET 心 肌 灌注 仿真 数据 

在 我 们 的 仿真 研究 中 ,我 们 选用 单 房 室 模型 ,如 图 
1 所 示 , 其 中 K(mLAmin/g),k(lmin) 为 2 个 标准 的 动力 
学 参数 。 左 边 红色 框 内 的 是 动脉 血 ,右边 蓝 色 框 内 的 是 
心肌 组 织 , 放 射 示 踪 物 以 K 的 速率 从 动脉 血 进入 心肌 
组 织 , 并 以 反馈 速率 k; 从 心肌 组 织 流 回 动脉 血 *"7。 


Uptake Rate (mL/min/g) 


Arterial Blood 


C(t) P 


Washout Rate (lmin) 
图 1 动态 ”Rb PET 心 肌 灌注 单 房 室 模 型 


Fig.1 One-tissue compartment model for kinetic modeling of 
“Rb MP PET imaging. 


这 两 个 速率 常数 与 动脉 和 心肌 组 织 中 的 放射 性 活 
度 成 线性 相关 。 动 态 过 程 可 以 通过 如 下 方程 描述 : 

ceo) = K,C,(t)- k2C (0) (1) 
FHC, OMC Qo BFE KAM AL LB PRB 的 浓度 。 
通过 测定 左 心室 (LV) 中 心 一 小 区 域 的 "Rb 的 浓度 获得 
Civ(0), 这 是 一 种 非 侵 入 方法 ,几乎 可 以 完全 恢复 动脉 输 
入 曲线 ,并 最 大 限度 减少 心肌 溢出 (C,(D)= Cry) ot 
且 , 如 果 定 义 五 为 组 织 中 的 血浆 体积 分 数 , 那 么 心肌 中 
的 放射 示 踪 物 浓度 可 以 表示 成 : 

C(t)=F,Ct) +0 —F Cn (t) (2) 
HP, FEO B 1 ZB SER, CRE NL AAR 
数 ,包括 : 

(a) 由 于 PET 采 样 时 的 部 分 容积 效应 从 血 池 进 入 心 
肌 的 

(b) 肌 肉 中 的 动脉 血 

对 方程 (1 ) 和 (2) 求 解 可 以 得 到 如 下 结果 : 

C()=F,C,()+(0-F) Ke" ® C,(t) (3) 
其 中 ,是 卷 积 算 子 。 


首先 为 多 个 区 域 (心肌 , 肺 等 ) 各 自 估 计 K,k; 常 数 
和 ,然后 求 其 平均 值 。 动 脉 血 放射 示 踪 物 浓度 Ci(? 通 
过 计算 动脉 血 采样 的 平均 值得 到 。 利 用 XCAT 体 模 *， 
根据 (3) 式 便 可 以 生成 一 套 动态 PET 心 肌 灌 注 图 像 。 如 
图 2 所 示 , 图 2A 给 出 了 K, 参 数 图 像 , 图 2B 为 采用 单 房 
室 模型 产生 的 时 间 活 度 曲 线 (TAC)。 我 们 在 2 min 内 采 
样 10 帧 (10x12 s), 这 是 心肌 灌注 PET 成 像 的 常规 采样 
协议 。 
1.2 PET 仿 真 断 层 成 像 

我 们 仿真 了 GE 的 RX PET/CT scanner 型 号 扫描 
器 呈 。 仿 真 器 产生 无 噪声 投影 数据 ,该 数据 经 过 放大 ， 
添加 泊 松 噪声 达到 真实 的 临床 噪声 水 平 。 仿 真 过 程 考 
虑 了 衰减 ,正则 化 的 影响 ,这些 影 响 都 包含 人 了 再 重建 
过 程 。 将 得 到 的 数据 分 别 采用 MLEM、 低 秩 征 罚 iit 
惩罚 和 L&S 四 种 方法 重建 。 其 中 MLEM 算 法 ,经 过 63 
次 迭代 得 到 (等 于 用 OSME 算 法 ,21 个 子 集 ,3 次 迭代 重 
建 吧 ) 的 动态 图 像 是 128x 128 x47 的 三 维和 矩阵 , 体 像 素 
大 小 为 3.27 mm x 3.27 mm x 3.27 mm. 
1.3 IKE EKA 

将 低 秩 矩阵 X 从 观测 矩阵 MN 恢 复 的 问题 ,如 (4) 式 
所 示 , 其 中 MM 是 包含 误差 的 矩阵 MM =X + E) 

min 5 // X-M //?+drank(X) (4) 


Herp rank(X {CASEMENT // + // NPL, EES 
数 ,用 来 平衡 第 一 项 ( 保 真 项 ) 和 第 二 项 (惩罚 项 )。 直 接 
求 (4) 式 中 目标 函数 的 极 小 值 是 很 困难 的 ,因为 里 面包 
含 秩 惩罚 。 为 了 简化 计算 ,将 秩 惩罚 用 核 范 数 包 X/ « 
表示 ， 

MX/ = do, (5) 
其 中 有 是 矩阵 X 的 奇异 值 。 经 过 这 样 的 松弛 处 理 , 低 秩 
矩阵 X 的 恢复 可 以 简化 为 ; 

min //X-M //?+2//X// + (6) 
1.4 L&S 模 型 重建 PET 动 态 灌注 图 像 

对 于 PET 动 态 心 肌 灌注 成 像 令 x 表 示 第 ! 帧 图 像 的 
示 踪 剂 活 度 ,1=1,2,…7,T 是 动态 图 像 总 帧 数 ,x 表示 第 
帧 第 个 体 素 的 示 踪 剂 活 度 ,j= 1,2,.…N ,入 为 体 素 的 总 
个 数 。 由 此 ,我 们 可 以 构造 矩阵 ; 


2 了 
XxX, X1 eee Xl 
x) x3 eee x 
2 2 2 
Mey E o E (7) 
1 2 T 
Xy XN eee XN 


其 中 ,XeR 代表 所 有 的 动态 帧 图 像 。 和 矩阵 XX 的 每 一 
列表 示 相 应 帧 图 像 的 示 踪 剂 活 度 , 每 一 行 表示 对 应 体 素 
AYES TANG EE HAAR (TAC). SEMEN SCI EAE SEER PA 
JERNE. FB Retail x a [ay fey BEA 
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图 2 KK 参数 图 像 和 采用 单 房 室 模型 产生 的 时 间 活 度 曲线 (TAC) 


Fig.2 K, parametric image(4) and TACs generated using one-tissue compartmental model(B). 


关 的 部 分 (例如 肌肉 组 织 和 肺 ) ,这 部 分 的 示 踪 剂 活 度 随 
时 间 变 化 缓慢 ,可 以 将 其 看 作 低 秩 部 分 (L)。 动 态 部 分 
代表 帧 或 区 域 之 间 变化 的 部 分 (例如 心肌 和 血 池 ), 这 部 
分 的 示 踪 剂 活 度 随时 间 快 速 变 化 ,是 稀 踊 的 (S)。 因 此 
我 们 在 动态 PET 心 肌 灌 注 成 像 上 同时 引入 低 秩 和 稀 玻 
惩罚 , 即 L&S 模型 。 求 解 过 程 可 以 表示 成 对 一 个 凸 优 
化 问题 的 求解 of。 

min /X-M //? +u AX) + uX) (8) 
其 中 ,pi A pe, ERRENA, AAEE ET 
项 ,将 其 用 p- 范 数 表示 ,定义 为 4(X)=(S Yo,")”, 用 来 惩 
罚 小 的 奇异 值 。 当 p < ] 时 ,该 惩罚 项 为 非 凸 的 范 数 , 例 
如 Lingala 选择 p=0.1 的 谱 惩罚 吓 。 这 里 ,我们 选择 
P=1, 将 其 简化 为 核 范 数 。g(X) 定 义 表示 时 间 空 间 的 全 
变 分 项 ,也 就 是 矩阵 在 x,y,t 方 向 的 梯度 的 1- 范 数 ,通过 
有 限 差分 估计 ,定义 如 下 : 

>| OY jor xe, (9) 
其 中 , =D, =D, AA AEE E 
AERE, D n Dp, DÆI x y A RAR 
1.5 最 优化 方法 

本 文采 用 split Bregman 法 "来 最 优化 (8) 式 中 的 
代价 函数 , 它 本 质 上 来 说 是 一 种 增 广 拉 格 朗 日 法 。(8) 式 
中 代价 函数 的 求解 可 以 通过 下 列 迭代 实现 : 


Xe = arg min||X-M-+2Z'||_ +p] x- E +r |+ 


B57 or xe -sivi 
EH! = arg min| X +L- E|} +2y/AilE|. GD 


S**!=arg min)? orx; -Si+ vil. 


C12) 
2 2.2|| ， 
+24s/B OS ds) | ,i=0,1,2 
1 
Latta x4 pe! (13) 
Vi avi +0 x8 i=0,1,2. (14) 


PLF aX GM (15) 


2 结果 
为 比较 前 面 提 到 的 不 同 重建 方法 的 效果 ,我 们 比较 
了 定量 评价 标准 。 均 方 误差 定义 为 : 


MSE Ror = Bias tor + Varror (16) 
其 中 ， 
ME Wa truey2 
(T= 
Biaskor = 一 - = 5 
yo (17) 
R N r yr \2 
Vans 1 2i E = 
ROI ! rue 
RS Gey (18) 


R 是 方差 实现 的 总 数 (我 们 的 仿真 中 R=30),N 是 
ROI 中 的 体 素 个 数 ,二 = 去 》)_ 叉 代 表 PET 图 像 在 第 j 
个 体 素 点 的 总 平均 值 ,x 代表 在 指定 ROI 第 j 个 体 素 经 
过 r 次 重建 ,x)“ 代 表 参 考 组 中 第 j 个 体 素 点 的 真实 的 
PET 图 像 值 。 为 了 比较 不 同 重建 算法 的 性 能 ,对 每 一 帧 
动态 图 像 计算 其 均 方 误差 (MSE), 如 图 3 所 示 , 横 坐标 
是 帧 , 纵 坐标 是 每 个 重建 算法 对 应 的 帧 的 MSE。 从 图 
中 可 以 明显 地 看 出 MLEM 方法 每 一 帧 图 像 的 MSE 值 
都 明显 大 于 其 它 3 种 算法 。 采 用 低 秩 惩罚 (L) 重 建 的 图 
RAY MSE 要 小 于 采用 稀 玻 惩罚 (S) 重 建 的 图像 的 
MSE。 我们 提出 的 L&S 方 法 重建 的 图 像 的 MSE 是 所 
有 方法 中 最 小 的 。 

为 了 提供 直接 视觉 比较 ,图 4 是 给 出 了 上 述 4 种 算 
法 重建 的 图 像 。 这 些 图 像 是 4D 动态 PET 心 肌 图 像 第 
26 层 断层 图 像 的 第 10 帧 。 从 图 中 可 以 清楚 地 看 到 
MLEM 算 法 与 稀 芍 惩罚 (S) 大 幅 增 加 了 品 声 水 平 , 低 秩 
惩罚 LL) 和 L&S 方 法 明显 降低 了 噪声 水 平 。 与 低 秩 惩 
罚 (LL) 重 建 相 比 ,L&S 重建 保留 了 更 多 图 像 特征 。 

此 外 ,为 了 评价 心肌 缺损 ,我 们 还 生成 了 左 心室 靶 
心 图 所 。 我 们 采用 美国 心脏 协会 的 17 段 模型 , 它 更 适 
合 评价 室 壁 运动 和 心肌 灌注 候 。 图 5 显示 了 不 同 重建 
算法 得 到 的 第 10 帧 图 像 的 17 段 靶 心 图 。 通 过 比较 ,我 
们 可 以 发 现 L&S 方 法 生成 的 靶 心 图 的 和 参考 组 的 真实 
靶 心 图 最 相近 。 我 们 可 以 清楚 地 看 到 与 其 它 重建 方法 
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图 3 不 同 重建 算法 对 应 帧 的 MSE 

Fig.3 Plots of MSE vs frame number for reconstructed 
images by MLEM (black), sparse penalty(green), Low 
rank penalty (blue) and L&S (red) algorithm. 


相 比 ,L&S 重 建 图 像 的 靶 心 图 无 论 是 在 正常 心肌 区 域 的 
还 是 缺损 部 分 ,与 参考 靶 心 图 的 偏差 都 最 小 。 
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3 讨论 

ASHE T PPS RS ASTI L&S) 
态 PET 心 肌 灌注 图 像 的 重建 方法 。 该 方法 可 以 表述 为 
一 个 凸 最 优化 问题 。 我 们 采用 split Bregman 算法 求解 
这 个 问题 。 我 们 用 单 房 室 模型 仿真 了 一 套 动态 PET 心 
肌 灌 注 图 像 ,利用 真实 的 动态 PET 心 肌 灌注 成 像 来 最 优 
化 正则 化 参数 。 本 文 提 出 的 L&S 重 建 算法 得 到 的 每 帧 
动态 图 像 的 MSE 是 前 面 提 到 的 算法 中 最 小 的 。L&S 重 
建 不 仅 降低 了 图 像 噪声 ,而且 保留 了 更 多 图 像 特征 。 此 
外 比较 不 同 算法 重建 的 图 像 的 靶 心 图 ,可 以 发 现 L&S 
方法 生成 的 车 心 图 的 灰 度 值 和 参考 组 的 真实 车 心 图 最 
相近 ,无 论 是 在 正常 心肌 区 域 的 还 是 缺损 部 分 ,与 参考 
靶 心 图 的 偏差 都 最 小 。 

总 之 ,我 们 提出 的 L&S 重建 算法 得 到 的 图 像 ,无 论 
是 在 直观 视觉 上 ,还 是 定量 分 析 上 都 优 于 另外 3 种 方法 。 


工 L&S 


图 4 从 左 到 右 依次 是 真实 图 像 (True) 以 及 通过 MLEM, RIEN (S) (RABE (L)ML&S 重建 得 到 的 图 像 
Fig.4 Raw image and reconstructed images by MLEM, sparse penalty, low rank penalty and L&S algorithm. 
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5 从 左 到 右 依 次 是 真实 图 像 (True) 的 靶 心 图 以 及 通过 MLEM , 稀 琉 惩罚 (S)、 低 秩 惩 罚 (L) 和 L&S 算 法 重建 得 到 


的 图 像 的 靶 心 图 


Fig.5 Polar map of the raw image and reconstructed images by MLEM, sparse penalty, Low rank penalty and L&S 


algorithm. All the polar maps were shown in the same color bar. 
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